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PREFACE. 

rjiHE lectures contained in this booklet were all delivered 
practically extempore. That they see the light of day, is 
chiefly due to the kindness of Mr. Gorakh Prasad, M.Sc., and 
Mr. Badrinath Prasad, M.Sc., Assistant Professors of Mathematics 
in the Benares Hindu University, of whom the former took down 
notes of the last three lectures and the latter of the first three. 
The copy for the press was prepared by me from these notes 
with the help of Mr. Mewa Ram Saksena, M.A., B.Sc,, of the 
Shyamsundar Memorial Intermediate College, Chandausi. Dr. 
Bibhutibhusan Datta, D.Sc., Acting Head of the Department of 
Applied Mathematics in the University College of Science, Cal- 
cutta, kindly helped me in correcting the proof sheets. I take 
this opportunity to express my thanks to all these friends and 
former pupils. 


Aprils 1924, 


CANESH PRASAD. 
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FIRST LECTURE. 

[Delivered on the 10th January, 1921.] 

INTRODUCTION. 

The chief reason for my choosing ^‘The place of partial 
differential equations in Mathematical Physios ” as the subject for 
these lectures is my wish to inspire in my audience a love for 
Mathematics. Before entering into details, however, I shall give a 
brief historical account of the application of Mathematics to natural 
phenomena. 

Mathematicians devote themselves to their subject chiefly with 
one of the two objects : either to study natural phenomena in con- 
nection with matter, space or time ; oi for purely aesthetic pleasure. 
A good illustration of the latter is afforded by their attempt to find 
all the prime numbers less than a certain given number. In ancient 
times people studied Mathematics mainly because of some mystery 



which they attached to it and not for its application in studying 
natural phenomena ; in this connection, we may note, for example, 
the doctrines of Pythagoras or the Sulvasutras of Apastamba. 
Among the earliest writers who devoted themselves to Mathematics, 
were Archytas of Tarentum (428-327 B.C.), Aristotle (384-322 B.C.) 
and Archimedes (287-212 B.C.). Archytas applied his knowledge 
of Mathematics in building up the science of Mechanics. After 
Archytas, Aristotle did something in this respect. But he had 
strange ideas about some elementary truths; for example, he 
supposed that for a falling body the rapidity of its motion depend- 
ed upon its weight, according to him a heavy body faUing more 
rapidly than a lighter one. Another erroneous idea which Aristotle 
had was that if a stone were thrown up at a certain angle, it 



2 DIFFERENTIAL EQUATIONS IN MATHEMATICAL PHYSICS. 

would go first upwards in a straight line at that angle and then 
fall vertically downwards, the path of the projectile consisting of 
two straight lines, as in the accompanying figure. Archimedes 
discovered the principle of the lever and the conditions of the 
floatation of a solid in a liquid. Galilei (1664 -1642) was the first 
to acquire clear ideas about the first two laws of motion. He also 
stated and proved that the path of a projectile was not a pair of 
straight lines but a parabola. After Galilei, Mechanics became 
gradually a science. A large number of results were discovered 
between 1642 and 1700. 

In the beginning of the 18th century a new epoch was created ; 
for, people began to think of natural phenomena other than 
ordinary motions, such as for instance the state of a hot body. 
This gave rise to the use of partial differential equations. Let there 
be a solid body which conducts heat. The amount of heat passing 


Y 



across a sheet at P perpendicular to the ^r-axis is required to be 
found out first. The amount of heat conducted will be propor- 
tional to the gradient of the temperature and is equal to 



which is the first germ of truth in the Analytical Theory of Hestt 
and was discovered about the year 1747; the partial differential 
equation governing the state of a hot body was discovered much 
later. Another important discovery was made in connection with 
the study of sound. People had already acquired some idea of 
harmony, and the formula of Taylor, 



VIZ. 



FIRST LECTURE, 




was roughly known to Aristotle and the ancient Rishis of India. 
In 1636 Mersenne gave the above formula as the result of his 
experiments. But all such results were at first more or less of an 
empirical nature. The equation arising in the problem of the 
vibrating string, viz. 

dt^ ~ P Sx^ 

was given by D’Alembert as late as in 1747 and was the first partial 
differential equation to be carefully studied. It was solved in 
T 

two ways. Putting — = a* , 

y = /(a; -I- at) + F{x - at ) , 

was the first solution given by Euler and D’ Alembert. Another, 
which was given by Bernoulli, is 


* . . TIttX 

y — ^ An sm cos 
1 t 


n-irat 

~ir * 


In appearance, Bernoulli’s solution is periodic whereas th(‘ 
former solution is non-periodic. Therefore the two solutions 
formed the subject of a most interesting controversy. Fourier 
showed that Bernoulli’s solution was also right and that, contrary 
to the views of Euler, it can be used for the case of an arbitrary 
displacement. 

It is the deliberate opinion of modern authorities that before 
the use of partial differential equations, there was no sound basis 
upon which the science of Physics could rest. Soon after the time 
when partial differential equations came to be applied to the study 
of natural phenomena, Fourier made the following statement 
which will remain true for ever : — 

“ The analytical equations, unknown to the ancient geometers» 
which Descartes was the first to introduce into the study of curves 
and surfaces, are not restricted to the properties of figures, and to 
those properties which are the object of rational Mechanics ; they 
extend to all general phenomena. There cannot be a language 
more universal and more simple, more free from errors and from 
obscurities, that is to say more worthy to express the invariable 
relations of natural things. 

“ Considered from this point of view, mathematical analysis is 
as extensive as nature itself ; it defines all perceptible relations, 
measures times, spaces, forces, temperatures ; this difficult science 
is formed slowly, but it preserves every principle which it has once 
acquired ; it grows and strengthens itself incessantly in the midst 
of the many variations and errors of the human mind. 

Its chief attribute is clearness ; it has no marks to express 
confused notions. It brings together phenomena the most diverse, 
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and discovers the hidden analogies which unite them. If matter 
escapes us, as that of air and light, by its extreme tenuity, if bodies 
are placed far from us in the immensity of space, if man wishes to 
know the aspect of the heavens at successive epochs separated by 
a great number of centuries, if the actions of gravity and heat are 
exerted in the interior of the earth at depths which will be always 
inaccessible, mathematical analysis can yet lay hold of the laws of 
the phenomena. It makes them present and measurable, and 
seems to be a faculty of the human mind destined to supplement 
the shortness of life and the imperfection of the senses ; and what 
is still more remarkable, it follows the same course in the study of 
all phenomena ; it interprets them by the same language, as if to 
attest the unity and simplicity of the plan of the Universe, and to 
make still more evident that unchangeable order which presides 
over all natural causes.” 

We have the following important partial differential equations 
that occur in the study of natural phenomena : — 


(1) 

d*F 

.a^v 




" do(? ' 



(2) 

dF^ 


B^V . 

d*F 

d*Fv 

dt " 

c* ^ 

da;* ■*“ 


Bz^) 

(3) 

B^V 

1 / 

d*F 

a‘F 

d*FN 

Bi^ “ 


da;* ^ 

ay^ * 

d2* J 

(4) 

B^V 


d*F 

- 0 . 


Bx^^ 


d2* 



B^V 


d*F 

■f X*F 


(6) 

Bx^ 

Bf 

^ d2»‘ 

= 0 ., 

(«) 

B^ V 

B^V 

d*F 



Bx^ 

Bi/ 

d2* 

= — 47r/u. 

/7\ 

B-V 


d^F 


\n 

^ Bt^ 

^ ^ at- 

' da;*’ 



As an illustration of the development of the rigorous method 
for the study of natural phenomena, I will sketch an outline of 
the history of the proof of Poisson’s equation, 

V»F=:~47r/>. 

The proof of Poisson’s equation, usually given in ordinary 
text-books, is not quite correct. The proof is based on enclosing 
the point (x , ;/, z) in a small sphere where the density p is taken to 
be constant. F=: F, -}- where F^ is the potential due to the small 
sphere and Fj is the potential due to the rest of the attracting 
body. 57*F|==0 and thus v*F= Now if you postulate that p 
is constant everywhere in the small sphere, then of course Poisson’s 
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proof can hold true. But in case p is variable, this proof will no 
longer hold. Green took up this question in 1828 and said that p 
is not constant but its variation is very small in the sphere and is 
therefore negligible. The first rigorous proof was published in 
1840 by Gauss who took p and its first partial differential ooeflS.- 
oients to be continuous. Now suppose F, the potential at P(a;, y, z), 
is due to a solid for which p at any point Q is 

/ ^ CO8^0 

logi 


where PQ = r, the radius vector and 6 is the angle made by the 
radius vector with the axis of x. In this case Poisson’s equation 
is not valid at P ; for. 


di/ 


00 , = — 20 and oo — oo — x 

or 


has no meaning. I invite your attention to Professor Petrini’s 
generalization of Poisson’s equation, viz. 


Lt 

h-0, 

kz=:0, 

2=0, 


r 1 s ^F(a;-f^,y,g) OFt 1 S z) dF > 

L M dx ^ ifc V dy % ] 


where A;, and I tend to zero without the ratio of any two of 
them tending to zero or infinity. There are various advantages of 
this generalization, but in my paper On the failure of Poisson’s 
equation and of Petrini’s generalization,” I indicated a case in 
which this generalization fails. Take for example 

P=coslogi; 

V 


then there is no limit in this case except for very special values of 
hy hy ly and the generalization as well as Poisson’s equation fails. 

The application of the modern Theory of Functions of a real 
variable makes the use of partial differential equations more logical. 
Also the application ot the Theory of Integral Equations is very 
useful. A typical integral equation is 

6 

(p{x)=f(x) 


ic(x, modi, 


I Bulletin of the Calcutta Mathematical Society, Vol. 8, 1918. 
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where </> is the unknown function to be determined. The first im- 
portant paper which laid the foundation of the Theory of Integral 
Equations as a systematic branch of Mathematics was published 
by Professor Ivar Fredholm in 1900 and was later published in the 
Acta Mathematica in an extended form in 1903. Integro-differential 
equations of a somewhat complicated type are useful in the “ Me- 
chanics of Heredity,” i.e. those branches of Mechanics in which the 
state at a particular point of time is dependent not only on the 
immediate past but also on the whole period antecedent to that 
point of time; as for example in the phenomena of Hysteresis. 
Integro-differential equations are equations which involve not only 
the unknown functions under signs of integration but also the 
unknown functions themselves and their derivates. A paper on 
the subject of integro-differential equations by Professor Vol terra 
appeared in the Acta Mathematical Vol. 35. The subject of integral 
equations is of essential importance in studying the phenomena of 



radiation. The law of Kirchoff, - = constant, cannot be proved 

without the help of integral equations. The first rigorous proof 
was given by Professor David Hilbert in 1913. The kinetic theory 
of gases cannot be rigorously studied without the help of integral 
equations (see Hilbert’s paper in the Mathematische Annalen, Vol. 
72, 1912). 

There can be no better instrument for studying natural 
phenomena than integral equations, partial and ordinary, and 
integro-differential equations. The truth of this remark can be 
easily seen by considering the problem of Dirichlet : To find a 
solution of inside a given region (Fig. 3) when, on the 

surface bounding the region, F=/ (x, y, z)^ a prescribed function. 
By the help of the integral equation 
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the quantity p can be found as the strength of a double layer whose 
potential is the required solution of the problem. Integral equa- 
tions are also most serviceable in studying the problem of the 
conduction of heat in a body. The well-known partial differential 
equation for the linear conduction of heat is 

da;*’ 

But this equation is invalid for some very simple cases and has to 
make way for a partial integral equation. (See G. Prasad’s memoir 
Constitution of matter and analytical theories of heat ” in the 
Abhandlungen der Kdniglichen Gesellschaft der Wissenschaften zu 
Gottingen, 1903). 
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[Delivered on the 11th January, 1921.] 


THE PARTIAL DIFFERENTIAL EQUATIONS OF VIBRATORY 

PHENOMENA. 

The most important partisbl differential equation met with in 
the treatment of ordinary vibratory phenomena is of the form 

a 

dfi \ dx* ^ V ^ SzV * 


Vibratory phenomena of electrodynamio origin are governed 
by the equations of electrodynamics : — 


4nW: 


dy dp 

dy dz* 



dy 


dx dy 


dt ^ dy dz* 
dt ^ dz dx * 


, ‘ dc dX 

dt dx dy * 

where the components of magnetic force at the point (x, y, z) are 
denoted by a. p, y, the components of magnetic induction by a, h, c, 
the components of electric force by X, F, Z, and the components 
of total current strength by u, t?, w. These six equations represent, 
in symbols, the laws of Faraday, and their solution reduces to the 
solution of a single partial differential equation of the type 


dt^ 


c» 


\dx^ dy^ ^ dz^ ) 


+ 4tncht . 


In the study of many kinds of elastic vibrations, we come 
across a partial differential equation of a similar type. In the case 
of the vibrations of an elastic plate, the motion is given by solving 
an equation of the form 
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a^w 

dt^ 




a* a* a^ 

where, as usual, stands for gp “*■ * 

I proceed to give you two examples to show how partial 
iifperential equations are utilized to determine the vibrations in 
3onorete oases : — 

First consider the motion of air in a rigid spherical envelope. 
Transforming the equation 

a^v 

a«» 

into polar co-ordinates, we have 


a<» 




1 ^ 

sin* 5 9^1* 


i 


Putting V=U (r, 0, <p) , cos kct, we find that the equation is 
satisfied by 

-v r 


where <p) is any surface harmonic of degree n. Certain bound- 
ary conditions are to be satisfied in every case. Here 

Jn*h(ka) 

' s/a "■ * 

a being the radius of the envelope. This gives the suitable values 
of k and thus we have an infinite number of solutions each satis- 
fying the differential equation and the condition at the boundary. 
The sum of such solutions 


J {k r) 

7=2 {9, COS kct 

».x Vr 

must also satisfy the differential equation and the boundary con- 
dition and gives the motion, the constants in /S„, being so chosen 
as to give the prescribed initial state when t is equated to zero. 

As a second example, consider the vibrations of a circular 
membrane. 

® Via** Syy 

is the equation of motion. On transforming it into polar co-or- 
dinates, we have 

a^v 


a^v 1 av 1 
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Proceeding as in the first example, 


BW 1 BU 1 B^U ^ 


The typical solution of this equation is 

U = {kr) (^4„ cos n $ Bn sin n d). 


Since at the boundary r=a, there is no displacement, U is 
zero there and therefore A; is a root of the equation 

Jn{ka)=i0. 

The motion is given by 

F = 2 cos ct (A^ ^ cos nB + ^ sin nB ) , 

the constants being determined by the prescribed initial state. 

Now I want to give you some criticisms, from the point of 
view of the theory of functions of real variables, on the usual 
method adopted for the study of vibratory phenomena. For the 
sake of simplicity and fixity of ideas, I will consider the simplest 
vibratory phenomenon, viz. the vibrations of a string. 

The partial differential equation of the motion is 

T 

where c* stands for — and may be taken to be unity if the units 

are suitably chosen. Let the string be fixed at its extremities 
a:=o, x:=ir and let it be set in motion by a displacement i/ = / (^). 
Then the motion at any subsequent time t is given by 

<ac 

y=- % i4„sin nx cos nt , 

»=sl 

the constants being chosen in such a way as to satisfy the initial 
condition . Therefore 

ao 

f{x) = 2 4„sin nx 


and 


An 



0 


fix) sin nx dx. 


The fiaws in the above method for determining the motion of 
the vibrating string were pointed out by me in a paper ^ “ On the 


I See Vol. 6 of the “ Bulletin of the Calcutta Mathematical Society, '' 1916. 
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ibratiixg string with an infinite number of edges.” The following 
re the two flaws : — (1) We always assume the existence of 


at* a? 


t every point. (2) It is taken for granted that if each term of an 
ifinite series satisfies the differential equation then the sum of the 
eries will also satisfy the equation. 

Both the assumptions are unjustifiable. For, a string may 
ave an edge so that at that point 

^ , what to speak of ^ , 

dll he non-existent ; also what is true of a finite number of terms 
lay not be true of an infinite number of terms. These objections 
proceed to elucidate at some length. 

Consider the string which is under constant tension T. Ac- 
ording to Newton’s second law, study the change of momentum 
f a small element of the string in the period of time 

(^, -T, f, +t). 

jet the oj-co-ordinates of the ends of the element be cc, - A, + 
?heh the change in its momentum is 


Xi t h 



® 1 —h 



dx . 


It must be balanced by the changf^ in the force which is 


ti+T 





)nly the assertion of this equality gives the true equation of the 
)henomenon as an integro-differential equation. 

What is usually done is to consider 2 ^ as getting smaller and 
mailer and also 2 t getting smaller and smaller. Then the appli- 
lation of Rolle’s theorem shows that the integro-differential equa- 
tion is equivalent to 


2h . 2t 



= 2t. 




x=x\*is%h 


vhere the ^’s are all proper fractions. 
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The ooaclusion, that 

^ OU**! t ' ' XssXi , 

is not always right ; if there is discontinuity then the above equa- 
tion may be invalid. 

Now I proceed to elucidate the second flaw by taking the 
initial state of the string to be y=zf{x), where 

f{x)^x forO^a;^^ , and / (a;) = ?r — .r for ^^x^t , 

The ordinary procedure gives the displacement y at any time as 

4 * ^ sin {2w— l)a; cos (2a- l)^ 

^ (2n-l)^ ' 

the assumption being that 

- - 1 ( -!)”■' sin (2w- l)a; cos (2w- 1)< 

Ot It 


and = 2 (- If’^ sin (2a~ l)a; cos (2a- 1)< 

OX*^ Tt wsbI 

and therefore the equation 


is satisfled. 
to which 


But the assumptions are unjustified because the series 


dt^ 


as well as 




is equated has no meaning, it being not convergent. Thus both 

^ 

de* ^ dx^ 


have to be investigated by a different procedure and their equality 
established by a different method. 

How to study then the phenomenon of the vibrating string ? 
This question was answered correctly by the late Dr. A. Harnack ^ 
in 1887. If 


y 




1 See his paper, “ Ueber die raifc Ecken behafbeten Sohwingungen gespannter 
Sait<en** {Mathemaiische Annalen, Bd. 29). 
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are finite and continuous then it was proved by Harnaok that the 
most general motion of the string is given by 

00 

y = 5 A„ sin nx cos nt> 

n=:l 

Proof : — 

Since is existent and finite, according to Fourier’s theorem 
ax 


y {x, t) must be equal to 


2 (t) sin nx, 

fiss] 


2 r 

I 

"J 


y (x, t) sin nx dx. 


[Here it may be noted that every continuous function is not 
expansible in Fourier’s series. This was first shown by Du Bois 
Reymond who gave an example of a continuous function the 
Fourier’s series corresponding to which is not convergent.] 
Integrating by parts 


sin nx dx 


-”1 


y sin nx dx^ — n* - o„ (0- 

A 


Since is continuous, we have 


y sin nx dx— \ ^ ^ » 


which equality is true because ~ is continuous. Therefore similarly 

at 


y sin nx dx- 1 ^ sin nx dx ; 


and this latter integral equals 


sin nx dx, 


because 




Hence 
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Therefore cos nt + sin nt. 

If the string is simply displaced initially, then from 

ao 

^ (An COS nt-^Bn sin nt) sin nx 
?» = 0 

we find that the B’s must be all zero. 

Therefore the general solution is 

00 

y- An cos nt sin nx. 


1 proceed to consider the case when there is an edge at a point 

d?/ 

P on the string so fchat at that point ^ is non-existent and conse- 

dx 

quently ^ is also non-existent there. 


Denote the ordinates by y^, in the portion APB and by 
y'p y\ in the portion AP'B. The usual phenomenon of wave mo- 
tion is that every wave propagates itself. Consider the small ele- 
ment of the string round P during a small period of time t. This 



element is ^ t where c is the abscissa of P ; the element undergoes 

(It 

a change 

^ dt'^ t dt s 

in momentum and the change is brought about by the impulse 

T 

idx dx S’ 

T 

acting for the period of time t. Therefore, since - = 1 , the units 
being suitably chosen, at P 

dt X dt dt ^ da? dx 


( 2 ) 
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As the string is continuous (c, 0=2/2 0 » therefore 

By\ ^ ^ + ^ 

dc * dt^ dt dc ' dt dt ' 


This gives 


gyi ^?yi \ 

rfA do? dx ) 


(3) 


In order that the two equations (2) and (3) be consistent, — 

dt 

must be equal to 1. Therefore the equation, which must be sa- 
tisfied at the edge in the place of the partial differential equation, 
which has no meaning there, is 



The question is now whether in the present case the general 
solution holds ; it will bo easily seen that the required condition 
at the edge is satisfied by the solution and therefore the answer 
to the question is in the affirmative. 

Cases of strings with three or four or any finite number of 
edges have been considered by Harnack. But the case of an in- 
finite number of edges was treated for the first time by me in the 
paper already referred to. As a simple example, it is proved there 
that 

^ G 

w = - S cos nt sin nx 

TT 

gives the motion of the string, vibrating with the infinite number 
of edges, the abscissae of which initially formed the set 


being equal to 




y 1 x 2 00 f . Wtt f 

according as n is even or odd. 


I will conclude this lecture by showing how the theory of 
integral equations may be applied in the study of the problem of 
the motion of a vibrating string. 

Let us consider the case of a string of variable density. Then 
in the equation 


p{^) 


at* a** 


(4) 


P is not constant. Lord Rayleigh could not integrate this equation 
except for very special types of p. 
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The general solution can be obtained by the help of an integral 
equation. For, put 

y = cos nt . U{x). 

Then the equation (4) becomes 


dx^ 




and it can be shown without difficulty that this last equation to- 
gether with the boundary conditions 

?7(0)=C/(ir) = 0 

is equivalent to the integral equation 

IT 

£7(ic) = y, U{() di, 

<1 

where G (i. x) = for 

TT 

Q x) = for $\x. 

TT 



THIRD LECTURE. 

[Delivered on the I2th. January, 1921.] 


THE PARTIAL DIFFERENTIAL EQUATIONS OF THE ANA- 
LYTICAL THEORY OF THE CONDUCTION OF HEAT. 

After the discovery of the experimental foundations of the 
elementary laws of conduction and radiation by Newton, Lambert 
and Biot, a good deal of time elapsed before the discovery by 
Fourier of the following differential equation as representing the 
phenomenon of the conduction of heat in an isotropic solid : — 

where is the coefficient of diffusivity. A number of special 
problems were solved by Fourier and his contemporary Poisson. 
For example, take the problem of the cooling of a solid sphere, 
whose surface is maintained at a constant temperature, say zero : 
the initial state being given by 

V(x, y, 2 , 0) =/(a;, y, z) . 

The particular case in which f(x, y, z) is a function of r alone 
was solved by Fourier; the general case was considered by 
Poisson. To determine the state inside the sphere at time we 
begin by transforming the partial differential equation into polar 
co-ordinates. Thus 


dt 


'-'7) 

< dr/ 

Putting 



F= < 



we have 



d 


1 d 

dr 


sin 6 do 




sin* 6 d<t>^ 

of which the typical solution is 

where (0, f) is a spherical surface harmonic of degree n. At the 
boundary r=a, V is zero ; therefore k must satisfy 

(ka)=0. 


(B, (p) 


2 
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Let the values of k he k^, k^, etc. Then the required solu- 
tion of the problem is given by 




the summation extending over all values of n from 0 upwards and 
over all values of m from 1 upwards and the constants in being 
so chosen as to satisfy the condition relating to the initial state. 

As another example, consider the state of heat of an ellipsoid 
of three unequal axes. Lam6 in 1837, Heine in 1839 and Mathieu 
in 1870 made important contributions to the advancement of the 
knowledge of the subject. Transforming the partial differential 
equation v®F=0 into ellipsoidal co-ordinates, the equation 
reduces to 

(“* !)■''*"-» (“'I:)' -» 

^vhere Re stands for \/(a‘^ + 0) (b^ + 0) (c^-^0). Taking 7 to be of 

ihe form B(\) E(/ii) E(v) where E{0) is a function to be determined, 
I’he equation will be satisfied if 



E(0)^A$ + B, 


vhere A and B are constants. The function E($) is called Lame’s 
'unction and the product ^(A) ^(/«) E(v) is called an ellipsoidal har- 
nonic. The stationary state of heat of the ellipsoid can be 
expressed in terms of a series of such products. As regards the 
ion-stationary state of heat, the partial differential equation of 
bhe phenomenon is 


iind this reduces to the form 




vhere 


*(»-.) (ft D’v}.” 






Mathieu made an unsuccessful attempt to solve the above 
equation. For the special case of an ellipsoid of revolution, Professor 
2). Niven’s results are of some interest. The most successful 
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attempt to solve the problem of the non-stationary state that has 
been made up to now is that of Mr. Bibhutibhuson Datta, who 
obtained, for the general case of the ellipsoid of three unequal 
axes as well as for the particular case of the ellipsoid of revolution, 
approximate results of very great interest and importance, one of 
the results being the correction of a mistake of Prof. Niven. (See 
the Bulletin of the Calcutta Mathematical Society, Vol. 8, and the 
American Journal of Mathematics, Vol. 41). 

I propose now to submit the use of partial differential equa- 
tions to a criticism from the point of view of one who wishes to 
be strictly rigorous. For the sake of simplicity and fixity of ideas, 
I will consider the conduction of heat in an infinite slab bounded 
by the faces a;= - tt, a;==7r, which are impermeable to heat, and I 
will take the initial temperature of the slab to be f(x) where f{T) 
is an even function of x. 

The partial differential equation of the phenomenon is 
dt ^ 

where stands for the coefficient of diffusivity and may be 
taken to be unity if the units are suitably chosen. Also because 
of the impermeability of the bounding faces, x=: tt, = we 
BV 

have -;^ = 0 f(^r = as well as for aj= - Thus the state of the 
Bx 

slab is given by 

V = + 2 e~ cos nx, tXO, — tt . 

1 

where 

TT 

= - I f{^) cos nx dx, 

0 

The above-mentioned treatment of the problem is the ortho- 
dox one and will be found in every text-book. In Fourier’s 
original papers the statement of the conditions is the same as 
above, with this difference that it is clearly mentioned that 

BV B^V 
Bt 

for every positive value of t including zero and for every value of 
X in the interval { - tt, n). 

In order to establish the defective character of the orthodox 
treatment, I take the simple case in which f(x) is given by the 
accompanying graph. 



20 DIFFERENTIAL EQUATIONS IN MATHEMATICAL PHYSICS. 


At O there is no differential coefficient and still by the ortho- 
dox method 

TT 4*cos(2»i— l)a; 

If we differentiate the above 
series term by term the series 
so obtained is uniformly con- 
vergent as long as / > 0 and 
hence we can differentiate in 
that manner for such values 
of t. Thus for f > 0, 

dF(0, * -{2n~\)H 

— — " '■ =: 2 ) 6 

dt 1 

Fig. 5. and, making t tend to 0, 

lim 3^(0, <) 

<=0 at ~ - 

Now according to the theory of functions of real variables, if 

lim SV{Q,t) 
t = 0 

is existent, it must be equal to 

lim F(0,<)-F(0,0 ) 

< = 0 ( 

(IV 

Therefore, at the time <=0 the value of at 0 is infinite 

at 

and also at that point, there being no differential coefficient of V 
with respect to non-existent ; consequently for a; = 0, # = 0, 

the equation 

BV d^V 
Bt 

has no meaning and the orthodox treatment fails. 

How to study then the phenomenon of the conduction of 
heat ? 

This question was ansAvered correctly for the first time in my 
memoir, “ Constitution of matter and analytical theories of heafc 
published by the Royal Society of Sciences of Gottingen in 1903. 
(See Ahhandlungen der Koniglichen QeselUchaft der W issenschaften 
zu Oottingen, Bd. 2, No. 1). In the theory, as remodelled by me, 
the equation of the phenomenon is shown to be not a partial 
differential equation but a partial integral equation, and attention 
is drawn to the simple principle that it is physical continuity 
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which is necessary and not the continuity of the form of an ana- 
lytical expression. T proceed to elucidate these features of the 
theory. 

First, as regards the equation of the phenomenon, this 
follows from the application of the principle of the conservation 
of energy in accordance with which the amount of heat passing 
into a volume must be equal to the amount of heat absorbed to 
raise the temperature of the mass of the volume. In the form 
of symbols, this condition may be written as : — 

t X 

(^) I 0) ^dx' , v 

- TT 


T{x, t) denoting the temperature. It can be easily shown that 
this integral equation is equivalent to the differential equation 


BT d^T 
Bt Bx^ 

BT B^T 

if — and are finite and continuous. For, by the mean value 
theorem, 

H 

J I ax^ 


Bxx 


I dt' 


and being proper fractions. Similarly 


aj=a;i-hei (a?a~aj,) 
+0a(fa~^i) 


x% 

J I T(x', i^)-T(x', <|) = 


»i 



aj=sa?l-hf»3 (a;a— *!)» 
t=ty +04(<a-«i). 


From {A), the two integrals given above are equal; consequently, 
proceeding to the limit by making ^2 - and - x^ tend to zero, 
we have 

BT^B'^T 
Bt ""doj* 


for a;=a;,, if these differential coefficients are finite and 

continuous. 

Secondly, T{x, 0) need not be derivable from the analytical 
expression, which equals T(x, t) when #> 0, by putting f =0 in that 
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expression. Thus, in the remodelled theory the convergence of 
the Fourier’s series 

00 

+ 2 a^^ cos nx 


has no place. On the other hand the continuity of T(x, t), i.e., 
the equality of 

lim V(x, t) to /(a;), 

<= +0 


follows at once from the continuity of f{x) ; and the only restric- 
tion imposed on f{x), in order that V{x, t) should be the solution, 
is the following : — 

(B) For each value of x, within the interval ( - tt, the upper 


limit of 


dx 


in the neighbourhood of ^ = 0is finite. 

BV . 


The necessity for the finiteness of ^ is obvious ; for, since 
BV . 

- A; ~ is the amount of heat which flows across an unit area of the 
Bx 


plane at x in unit time, if this amount becomes infinite it will 
be absurd. Now for ^ > 0, the series for V (a?, t) being uniformly 
convergent and the series obtained by differentiating it being also 
BV 

uniformly convergent, is finite. Hence we have to consider only 

the neighbourhood of <=0. 

The equation {A) is satisfied by 

T(x^ t)—V(x, 0 , ^ > 0 

as long as f(x) is finite and integrable. For, for ^>0, the series 
dV 

for and V can be integrated term by terra, the series concern- 


ed being all uniformly convergent. 

BV 

As regards (H), it can be satisfied even if makes aninfinit< 
number of oscillations as t tends to 0 ; for example, if 

X 

0 

av , 

makes an infinite number of oscillations with indefinitely 



increasing amplitudes as t tends to zero. Thus theoretically the 
slab may, at certain points, become hotter and cooler an infinite 
number of times. 
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I conclude this lecture by pointing out that a number of in- 
teresting problems in the theory of the conduction of heat can be 
best solved by the use of the theory of linear integral equations. 
For example, one of the earliest applications of that theory is by 
Professor Picard, to the problem of the cooling of a ring. The 
partial differential equation of the phenomenon is 


— 


c® 


d^V 
* da;* 


>-6*F. 


This differential equation can be solved in terms of the solu- 
tions of an ordinary differential equation of the second order which , 
in turn, is equivalent to an integral equation of the type 

U(i) = F{i) ^ ^ ’ 

where U is the unknown function. 



FOURTH LECTURE. 

[Delivered on the 3rd February, 1931.] 

THE PARTIAL DIFFERENTIAL EQUATIONS OF THE THEORY 
OF GRAVITATIONAL ATTRACTIONS AND THOSE OF 
ELECTROSTATICS AND MAGNETOSTATICS. 

In the three lectures which I delivered last month, T first gave 
a historical review of the present theory of the partial differential 
equations of Mathematical Physics together with some criticisms 
and then dealt in succession with the equations of vibratory 
phenomena and of the analytical theory of the conduction of heat. 

I wish to preface the three lectures which I intend to deliver 
in this month by explaining to you the two kinds of Mathematical 
Physics. This cannot be better done than in the following words 
addressed by Professor Vito Vol terra of Rome at the inauguration 
of the Rice Institute of America in 1915 : — ‘‘ There are tw^o kinds 
of Mathematical Physics. Through ancient habit wo regard them 
as belonging to a single branch and generally teach them in the 
same courses, but their natures are quite different. In most cases 
the people who are greatly interested in one despise somewhat the 
other. The first kind consists in a difficult and subtle analysis 
connected with physical questions. Its scope is to solve in a 
complete and exact manner the problems which it presents to us. 
It endeavours also to demonstrate by rigorous methods statements 
which are fundamental, from mathematical and logical points of 
view. 

I believe I do not err when I say that many physicists look 
upon this mathematical flora as a collection of parasitic plants 
grown to the great tree of natural philosophy. But this disdain 
is not justified. In the evolution of Mathematical Physics these 

researches are to play an ever increasing part 

“ The other kind of Mathematical Physics has a less analytical 
character, but forms a subject insuperable from any consideration 
of phenomena. We could expect no progress in their study 
without the aid which this brings to them. Could any one imagine 
the electro-magnetic theory of light, the experiments of Hertz and 
wireless telegraphy, without the mathematical analysis of Maxwell 
which was responsible for their birth ? ” 

To-day I shall confine myself to the first kind of Mathematical 
Physios and the lecture will be mainly destructive in character. 
But to-morrow I shall speak on the constructive side of things, 
when I shall take up the second kind of Mathematical Physics. 
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The following are the most important equations of gravita- 
tional attractions, Electrostatics and Magnetostatics with the 
respective years of their discovery : — 


,,, d^V d*F d»F ^ 

^ 


(f being the surface density. 
a»F d»F a*F 

p being the volume density. 
(4) IF^o“W^-o--47rcr, 


(Laplace’s equation, 1781) 

(Equation of Coulomb and 
Poisson, 1811) 

(Poisson’s equation, 1813) 
(Helmholtz’s equation, 1853) 


(f being the strength. 

[The exact form of the equation (2) was given by Poisson in 
1811. The last equation gives the potential difference at the 
surface of a double layer of strength tr.] 

Let us take up Poisson’s equation 

S'ly S^V S^V 
dx^ ^ dy^ ^ dz^ 

The proof as given in text-books postulates p to be a constant 
throughout a sufficiently small sphere round (x, z). Let u and 

U denote respectively the potential due to the small sphere 
round {x, y, z) and that due to the remaining mass. Then 


V^U^u, 


By Laplace’s equation v ^ (7 = 0 
and consequently 

Taking p to be constant 

w = 27r/> (a^-ir*) 

at a point (^, v> 0 in the sphere distant r from (x, y, z),a being the 
radius of the sphere. Differentiating two times with respect to 
’/j £ we have at the centre 


r»M=-?^(2+2 + 2)=-4»rp. 

O 

In 1828, Green gave a proof more or less like the above. The 
first really rigorous proof was given by Gauss in 1840. Gauss’s 
proof, which T will reproduce presently, is based on the supposition 
that p is finite and continuous together with its three partial 
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differential coefficients of the first order. In 1882, Professor Holder 
succeeded in establishing Poisson’s equation on the supposition 
that 

I V. O-pion, y, *) I 

where m and A are positive constants and r is the distance between 
(f, Vi 0 and {x, y, z). The next investigator to impose still less 
restrictions on p was G. Morera who proved in 1887 that, if 



is finite, then Poisson’s equation is valid. In 1891, Kronecker 
attempted to establish the equation with lesser restrictions than 
the preceding investigators. From 1891 to 1900 practically nothing 
was done to improve the proof of Poisson’s equation. In 1900 and 
the next ten years. Dr. Henrik Petrini tried to make the proof as 
little restrictive as possible. Some of my papers deal with the 
criticism of Petrini’s work. Petrini ’s chief result is that if 






pdr 

r 




exists and is finite, then ^ exists and is finite, p standing for 

cos 9 where 9 is the angle between the radius vector and the axis 

of X ; with similar conditions for and 5-7 ; the additional con- 

oy oz 

dition for the validity of Poisson’s equation being the continuity 
of which condition, however, has not been sufficiently emphasized 
by Petrini. But cases may occur in which p is discontinuous. 
For example, I pointed out in a paper published in the Philo- 
sophical Magazine in 1917 that, in the case in which p is 


cos ^ log , Poisson’s equation is invalid. Instead of cos ^ log ^ ) 

we may take, without making any material difference, cos ^ log - ^ 
plus any constant or differentiable function. If, then, you take a 
solid sphere of density 2 -f cos ^ log , say, and suppose it divided 
into thin concentric shells, then, by Newton’s theorem, the attrac- 
tion due to a shell is zero at a point inside it and therefore 
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;r" = 0 at the centre. The mass of the shells between the centre and 
dr 

the point (f, 17, f) is 



47r 


^{2 + cos (log i)} 


dt. 


0 

Since the attraction at (f, v, t) is the same as if the whole of 
this mass were concentrated at the centre, the attraction is 


r 

j + (log J)} 


and is along the radius vector. Now at the centre 


du 

— = hm 


dr^ 


r = 0 


r 

= ^ I 2 + cos log i I 


r=0' 

47r r 


"r L“ [ - T - 5 ) }] ■ 

As r is indefinitely diminished, - tends to <» and cos ^ log ^ ^ 

d^u 

makes an infinite number of oscillations. Therefore VT is non- 
existent at the centre. Moreover the radius vector may be taken 
in any direction. Hence, at the centre, 

d^u d^u d^u 
dx^ ’ dy^ ’ dz^ 


are all non-existent and consequently the equation of Poisson is 
meaningless. 
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For the sake of its historical importance I proceed to give 
you Gauss’s proof of Poisson’s equation : — 

Let = Then the ic-component of the attraction at 

the point {x, y, z) is 


X(x, y, «) = 




{(!-»)» + yY + (f- 


( 1 ) 


where dt is the element of volume of the sphere and the integral is 
taken over the whole sphere. 

Also 


X{x ^ e, y, z) = 




{i-x-e) f(i, rj, g) dt 
{($-x-eY + {r,-yY + {C-zY}^ 


( 2 ) 


the integral being taken throughout T, the original volume occupied 
by the sphere. Suppose the sphere to be displaced a little to the 
left so that its centre is at i. e. ( - e, 0, 0). Let denote the 
volume occupied by the sphere in its new position. Then, as is 
clear from the accompanying figure (7), 



Fig, 7. 




Now, put 
Then 




X(a; + e, y, «) = 






(C-x)M + e, v,0dt 
{ii'-xf + (y-yY + (i-zf}^ 


(3) 
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Also 


y,z)-X{x, y.^) ^L|(UA)-(^|+X|)]- = 


e 


A, 


where I is the integral in (1) taken over x the same integral 
taken over d, the integral in (3) taken over and Ai the same 
integral taken over 

The value of an element of volume in 6 and that of an 
element in are respectively e cos a dS and -- e cos a dS where a is 
the angle made by the external normal with the x - axis. 
Therefore in the limit, 


L^={(i-X) f{(, v,0-~s^ dS 

taken throughout the surface S of 
also in the limit 







dt 


taken throughout T^, . 


Therefore 

f Izf ^ ^ i 

dx dx^ ) r® ) r* 

h 6’ 

and similarly for 

d^u - 
— and :rT- • 


pdSi 


Therefore 


Since 



To S 

dt = r* sin BdOd<pdr, 


the volume integral in the above equality will tend to zero as the 
radius of the sphere is made to tend to zero; and similarly 


J ^ will approach -4ir/ (a;, y, z). 


Thus the equation 

V * F = — 47rp 


holds. 

I proceed now to consider the generalization of Poisson’s 
equation given by Dr. Petrini. Dr. Petrini was the first to note 
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clearly that the continuity of p was not sufficient for the validity of 

cos ^ 0 

Poisson’s equation. Taking p= he showed ^ that Poisson’s 

log 7 

d^V d^V 

equation fails because ^ is + <» while and ^ are both equal 
to - 00 . Dr. Petrini was thus led to the following generalization : — 


lim r 1 f dV (a; 4-^1, y, z) dV (x, y, z) ( 1 f dF z) 

L ^ ^ ^ ^ ^ ^ % 

fiV(x,y,z)\ 1 fdV(x,y.z-^h^) OF («, t/, 2 ) ^ 
Jz 


provided that while A,, are made to tend to zero in any 
h h h 

manner, , r* , ir should tend neither to zero nor to infinity. Dr. 

tl^ /tg 


Petrini proved the validity of his generalization for p = 


cos* 6 

log 7 


and, 


indeed, for every case where p is continuous. But he did not con- 
sider the oases in which p has a discontinuity of the second kind. 
Such cases have been studied by me in the paper referred to above. 
I have proved, for example, that when 


1 ^ 
r dr 


p = OOS ^log , 

(log tan" y) 


which shows that the limit in the generalization of Petrini will 
exist and be equal to the limit of 


-47r 

^/lO 'I 


J 


COS 


^ log ^ + tan + cos ^log ~ + tan ~ ^ 


cos 



-1 

+ tan 



if the latter limit exists. Now put and Then 

1 

denoting log + tan g by log the expression within the crook- 
ed brackets becomes 


I For another method than Petrini’s, see my paper in the PhiL Mctg^. Vol. 
. 34 , 1917 . 
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COS (log H) + cos ^log H + log + cos ^log H + log , 
i. e., cos (log i^) I 1 + cos ^log + cos ^log | 

- sin (log H) I sin ^log + sin ^log | . 

Now the only way in which this expression can have a limit 
is by the vanishing of both 

{l+cos (logl) + oos and I sin ^lor; + sin (-I)S 

which is possible only when a and p have respectively the forms 

g2m^±2»/3 g2n^q:2,r/3 

m and n being integers. Thus, it is proved that the generaliza- 
tion of Petrini does not hold for every discontinuous distribution. 

Another generalization of Poisson’s equation was given by 
Professor C. W. Oseen of Upsala. As I have already taken a 
good deal of your time this evening, I will take up the consider- 
ation of that generalization to-morrow. I will conclude the 
lecture by discussing briefly the equation of Coulomb and Poisson 
and Helmholtz’s equation. 

I have studied the equation of* Coulomb and Poisson in three 
papers one of which appeared in the Rendiconti del Circolo Mate- 
matico di Palermo in 1917, another in the Philosophical Magazine 
in 1918 and the third in the Bulletin of the Calcutta Mathematical 
Society in 1919. If round P, any point on the surface, a small disc 
be described lying on the surface, then for investigating the validity 
of the equation 

- — 4. _ = — 47r(r , 

we need only take into consideration the part of F due to the disc. 

Among other results, I have proved that, denoting by r the 
distance of any point of the disc from P, if 


and, further, if 


(T = cos x(r) when x(r) 


X {r) log 


dV dV 

then ^ , — are both non-existent and the equation of Coulomb 


fln, ’ 


and Poisson fails. 

Results similar to those mentioned above have been obtained 
by me for the validity of Helmholtz’s equation in a paper which 
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appeared in the Proceedings of the Benares Mathematical Society 
for 1920. I have shown, for example, that, if the strength 

<r=coalogi 

r 

then TT + o are both non-existent and therefore the equation 

fails. 
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[Delivered on the 4th February, 1921.] 

THE PARTIAL DIFFERENTIA!. EQUATIONS OF THE 
THEORY OF GRAVITATIONAL ATTRACTIONS, AND 
THOSE OF ELECTROSTATICS AND MAGNETO- 
STATICS {Continued). 

Yesterday I promised to lay before you the generalization of 
Poisson’s equation given by Professor C. W. Oseen. The general- 
ization ' is fairly simple and will prove as useful as Poisson’s 
equation. Tt is this : — 



where dv is an element of the volume O 5 , ds an element of the 

dV 

surface 8 which bounds the volume Qf; and ^ is the differential 

coefficient with respect to the normal taken outwards. The only 
condition imposed here on p is that it should be finite and 
integrable. 

It will be remembered that in proving the theorem 

V^F=:-47rp 

we have to impose more restrictions on p than that it should be 
finite and integrable ; for example, in Gauss’s proof, it is supposed 
that p is continuous and also that p^,, p^, p, exist and are conti- 
nuous. We see at once how much more simple Professor Oseen ’s 
generalization is. 

Starting from the equation, known in Integral Calculus as 
Green’s theorem, 

one can deduce the modified theorem of Green, viz., 



i See Oseen’s paper, Ueber einen Satz von Green und.ueber die Definitionen 
von rot und div ” (Rendiconti del Oircolo matematico di Palermo^ t. 38, 1913). 

3 
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where stands for p{P') and V' for F(P'), P' being the variable 
point (^, 7/, c), P the fiyed point (x, y, z), and r the distance PP'. 
This result follows from Green’s theorem which holds when 
are existent and integrable. To show this, I take U to 

pp-, • Now Yp becomes infinite at the point P. In order 

to make Green’s theorem applicable, we enclose the point P within 
a small sphere il' of radius a. Inside the region bounded by this 
sphere and S, U and its differential coefficients are finite and 

continuous. Hence, since ~ satisfies Laplace’s equation, Green’s 
theorem applied to this region gives 



where ds' is an element of the surface S' of the sphere O'. Now 
let a tend to zero, then, dv being of the order r*dr, the second 
term in the left-hand side of the above equation is negligible. 
Similarly, as ds' is of the same order as r^, the third term in the 
right-hand side is negligible. Also v^V' - according to 

Poisson’s equation. Hence we have in the place of (1), 



which becomes, in the limit when a tends to zero, 



This result is the modified theorem of Green and it says that 
the value of F at any point is the potential due to (1) a distribu- 
tion of volume density p', (2) a certain distribution on the bound* 
ing surface and (3) a certain double layer on that surface. Thus 
we see that if there is a solid with a known volume distribution 

and prescribed values of F and ^ on the bounding surface, then 

by applying the modified theorem of Green we can at once put 
down its potential F at any point. 

The question arises : Can we deduce the modified theorem of 
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Oreen by using Oseen’s generalization instead of Poisson’s equa- 
tion % The question has been answered in the affirmative by Oseon 
who has deduced the modified theorem of Green by applying 
vector analysis and by employing the following theorem : — 

If <f> and are two functions, finite and continuous together 
with their first differential coefficients, then 


I (grad (j> X grad 4 /) dv 




The procedure adopted by Oseen is rather complicated, sinc(.* 
the second differential coefficients of (l> may be non-existont, but 
the deduction is an epoch-making result. 

The equation of Laplace, which 1 have frequently mentioned 
in some of my previous lectures is intimately associated with 
problems in Electrostatics and Magnetostatics. Thus we come 
across this equation when considering the case of a solenoidal 
distribution of electricity and the most important problem associa- 
ted with the equation is Dirichlet’s problem, viz. to find the 
potential at points inside a region when the value of the potential 
is given on the boundary. 

With reference to Dirichlet’s problem, two theorems have 
been formulated and proved under certain general conditions. 
These theorems are : — 


( 1 ) (The Theorem of Uniqueness) : The solution of Dirichlet’s 
problem is unique. 

( 2 ) (The Existence Theorem) : The prescribed value on the 
boundary is arbitrary and corresponding to each value there exists 
a solution. 

That both the theorems hold, was at first guessed from certain 
physical considerations. For example, in the case of a conductor 
connected with the earth and influenced by external electric 
masses, the value of the potential of the induced distribution on 
the surface is prescribed, being equal in magnitude, and opposite 
in sign, to that due to the masses, and the uniqueness and 
existence of the potential is obvious to a physicist whatever the 
masses may be. But attempts were made to base the two 
theorems on something more rigorous. 

Dirichlet formulated the principle, known as Birichlet’s 
principle, that, of all the functions, F, which inside the region 
bounded by a closed surface remain continuous together with 
their first differential coefficients and take a prescribed value on 
the surface, there must exist one (or more) for which the integral 
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taken over the whole region has a minimum. Accepting this 
principle, we get from the calculus of variations the result that V 
must satisfy Laplace’s equation. 

But Weierstrass pointed out in 1860 that Birichlet’s principle 
was not always valid and that all that could be said was that the 
integral has a lower limit and that it may not have a minimum. 
In 1899, Hilbert established the validity of Dirichlet’s principle 
on the supposition that the prescribed boundary value is con^ 
tinuous together with the differential coefficients of the first two 
orders. 

Assuming that in a given case the solution of Dirichlet’s 
problem exists and is unique, the question arises of actually 
obtaining this solution, I propose to devote the remaining time 
to-day to the discussion of a number of methods. These are (1) 
the use of Green’s function, (2) the use of fundamental functions 
based on a specially suited co-ordinate system, (3) the mithode dii 
halayage, (4) the method of arithmetic means of Poincar^ and (5) 
the use of integral equations. 

Adopting the procedure used in deducing the modified theorem 
of Green from Green’s theorem but taking U to be ^ instead 

of i , we obtain for V(P) the simple expression 





here O is regular everywhere inside the region, where it satisfies 
Laplace’s equation, and equals - ^ at every point on the surface ; 

O -4- ^ is called Green’s function for the particular region. The 

problem of obtaining the solution of Dirichlet’s problem has been 
thus reduced to that of obtaining Green’s function which may 
be interpreted as the potential due to a unit of electricity at P 
plus the potential due to the distribution induced on S by this 
unit when S is connected with the earth. In some cases Green’s 
function is easily obtained by applying Lord Kelvin’s theory of 
images. For example, let ^ be a spherical surface of radius 
a. Then the value of Green’s function at any point Q is 

g 1 1 

RP'Q*PQ‘ 

where R denotes the distance of P from the centre of the sphere. 
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P' the inverse of P with respect to the sphere ; in fact 5 units at 

P' form the image of a unit at P. 

I proceed now to consider at some length the second method, 
viz. that in which use is made of fundamental functions based on 
a specially selected system of co-ordinates suitable for the parti- 
cular boundary surface. 

As an example of this method, consider the case in which the 
bounding surface is a sphere. Then, transforming v*F = 0 into 
the polar co-ordinates r, B, r=a being the bounding spherical 
surface, we have 


cir \ dr f ^ sin 0 dO dO/ sm*6^ 


We start by obtaining a typical solution of this equation in 
the form of a product of three factors, one a function of r alone, 
another a function of $ alone and the third a function of </> alone. 
Such a solution which is the desired fundamental function, is easily 
found to be 


where 


cos mp 
(cos B) or 

sin mp 


» 1 .2x2 

2"n. ! j " ^ 

dfA 




1)-. 


The solution of Dirichlet’s problem is given by 


F= 2 2 coamf -t- sinm^), 

?) = (! »j = 0 

where the and P’s are constants which have to be determined 
by the condition that on the boundary, r=a, F takes the prescrib- 
ed value, say f{ 6 , p). Now in actual practice the determination of 
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these constants is not an easy affair. Gauss had given * a method 
for expanding any arbitrary polynomial (in the co-ordinates of a 
point on the spherical surface) in a series of spherical harmonics. 
If <}>{x, ijy z) be the polynomial of degree n in x, y, z where 

then Gauss put 

1>(x, y, 8n-.i + , 


8^ standing for a solid spherical harmonic of degree m, and, by 
operating on the above equation with v*, etc., obtained a 
number of simultaneous equations by solving which the ^’s were 
obtained. But no com})act and general formula for writing down 
the expansion had been obtained until about twenty* years ago I 
gave it in a paper published in the Messenger of Mathematics, Vol. 
30, in 1900. 

The result is 


y. ^) = (- 

where 


1 ) 


r 2n-2Mi I r 2 u ^ l~i 

(•2n~4;.+l) r |r, ^6(a, y) | ^ 

2/^ {2/u - 2) . . . 2 . (2a -2/1 + 1) (m - 2/x - 1) . . .3. 1 




/3 = £ , 
^ % 




d 

^ 2 ’ 


\ a — 




da 


The above-mentioned result was obtained by me at first by a 
method different from that of Gauss whose method was applied by 
me to get the result in a paper published in the Mathematische 
Annalen, Vol. 75, in 1912. 

If / {0, <p) is expressed as a function of x, y, z, the co-ordinates 
of a point on the boundary, that function can, in its turn, be 
expressed as a sum of polynomials each of which may be expanded 
by my formula, and thus an expansion for f[0, <t>) can be obtained. 
With the help of this expansion, the solution of Dirichlet’s problem 
can be put down at once as was done by me in my paper of 1900. 

As a second example of the method of fundamental functions, 
consider the case in which the boundary is a cylinder. Then, 
transforming v*F = 0 into cylindrical co-ordinates ik 2 :, p — a 
being the boundary, we get 

d^V 1^ l_d^ ^!Z = 0 

V ^ * p^ “ ■ 


The typical solution of the above is 

±kz GO& mp 

e ikp) or 

sin 


I See his “ Gesammelte Werke**, Bd. 5, p. 360. 
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this is the desired fundamental function, and the solution of 
Dirichlet’s problem is given by 

F = QOS m<f> + sin mi>^ e + D„,m e jjm , 

where the constants, the ^’s, B’s, and D’s are determined by 
the boundary condition. 

As the third example of the method of fundamental functions 
consider the case in which the boundary is an ellipsoid. Then, 
transforming v^F = 0 into the ellipsoidal co-ordinates A, /i, r, A = 0 
being the given ellipsoidal boundary 


we get 


where 


Re = ^ . 


We suppose, as usual, that a typical solution of the above 
ei^uation can be obtained as a product of three factors, each being 
a function of one co-ordinate alone. Thus we get as the funda- 
mental function 

E{\) E{n) E{v), 

where 

E{e)=z(mO*7i)E(6) , 

which is Lame’s equation and can be solved in series The 
required solution of Dirichlet’s problem is given by 

Em,,t (A) (/•) (r) , 

the A’s being constants which are determined by the condition 
that, for A = 0 , V has the prescribed value, say F(/s >'). Many 
attempts have been made to obtain a compact and general formula 
for expanding » ) in ellipsoidal harmonics similar to that given 
by me for expansion in spherical harmonics, but so far the attempts 
have proved to be unsuccessful. 

The question of the determination of the most general system 
of three orthogonal surfaces, namely, that of the confocal cy elides, 
for which the fundamental function is the product of three factors, 
as in the cases considered above, was investigated by Professor 
Klein and his pupil, the late Professor Maxime Bocher. (See 
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Bocher’s book, “ Reiheneutwicklungen der Potentialtheorie 
Leipzig, 1894.) For the case of an arbitrary bounding surface, 
Poincare established in 1896 the existence of fundamental func- 
tions, which, however, are not necessarily expressible as the product 
of three factors, each factor being a function of only one co- 
ordinate. (See his paper in the A^a Mathematical Vol. 20.) 

Unlike the previous method, the meihode du halayage of 
Poincar^ ‘ uses approximation functions each of which has the 
prescribed value on the bounding surface and each satisfies 
Laplace’s equation more approximately than the one preceding it. 
Briefly, the method may be outlined as follows : — 

Let if be a sphere which encloses the surface S and the region 
ilf; bounded by it and let the prescribed value f(x, y, z) be such 
that / together with its differential coefficients of the first two 
orders is continuous in K ; also assume that in K, and let ({> 

denote 


A' 




. v y .dt . 


Then obviously v ‘^(|)~ /) = 0 in . Now construct an enu- 
merable aggregate of spheres if^, iTp K.^, . . . such that each sphere 
lies in Os and fills it in such a manner that each point of Og is 
inside at least one sphere and the spheres become indefinitely 
large in number only in the neighbourhood of S, Then, if the 
sequence of functions f/,, ... is determined by the condi- 
tions that (1) = ^ fA) equals </> in K - (2) 

^*f7, = 0 in if, and U, equals r/^iniC ~ (3) = 0 in and 

equals in etc., A U, A . . . and the sequence 

Up . . . converges uniformly to a regular potential-function 
in Qs : the required solution of Dirichlet’s problem is 


n = 00 




The case in which y V >0 can be dealt with by making it 
dependent on the case considered above by suitable superposition. 

The method of arithmetic means as given by Neumann * 
postulates the convexity of S and determines the solution of 
Dirichlet’s problem as the potential of a double layer. If 

S S 


* See his paper, “ Sur les equations aux d6riv6es partielles de la Physique 
Math^matique {American Journal of Mathematics^ Vol. 12, 1890). 

2 See his booklet, ** Untersuchungen uber das logarithmische and New- 
tonsche Potential,'* Leipzig, 1877. 
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8 


then it has been proved by Neumann that the series 

- }l:0h 

s 

converges in Us , satisfies Laplace’s equation and, with the approach 
of the point P, at which the potential is taken, to the boundary, 
tends to /. Here C is a constant and equals the value of 
lim on the surface. 

n=oo 

The method in which the theory of integral equations is used 
is allied to the above method inasmuch as the solution of l)>richlet’s 
problem is obtained as the potential due to a double layer. Let 
the strength of the layer be />, then it is known that 

We -"W i= — 4:irp(P), and 

8 

where We or is the limit to which the potential W of the 
double layer tends according as the point P on the surface ivS 
approached from the outside or from the inside, P, M, 'P having the 
same meanings as in the figure on p. 6. Adding the two equations, 
we have 

2W, = - (P) + 2 j p(M) ^ d, . 

s 

But We is equal to the prescribed value /(P). Therefore p is 
given by the integral equation 

s 
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[Delivered on the 5ih February, 1921.] 

THK PARTIAL DIFFERENTIAL EQUATIONS OF HYDRO- 
DYNAMICS AND THOSE OF ELECTRODYNAMICS 
AND THE THEORY OF ELECTRONS. 


I propose to speak to-day first on the partial differential 
equations of Hydrodynamics. Less than three hundred years 
ago, Hydrodynamics was not a science at all. In fact the first 
theorem ever discovered in Hydrodynamics is what is known as 
Torricelli’s theorem. This was first stated in the year 1643. 
Torricelli was one of the greatest engineers of his day. He found 
from observation that the velocity of water issuing from a small 
hole in a vessel of water is the same as that of a particle falling 
from the level of the water in the vessel, i.e., 2gz, where q is the 
velocity of water flowing out and z is the depth of the hole below 
the water level. The first attempt at a mathematical proof of 
this, was made by Newton. Ho pointed out that two important 
points were to be kept in mind. ( 1) The fluid particles in any 
horizontal section of the jet will always remain in a horizontal 
section. (2) As the section increases in size, the velocity in the 
section becomes less in the inverse ratio. This was a crude form of 
the equation of continuity. From considerations like these Newton 
satisfied himself that the theorem of Torricelli was correct. He 
published his arguments in the Princi'pia in 1687. But the 
theorem was first proved completely by D. Bernoulli in 1738. 
In 1755 Euler gave the equations of motion 

Du ^ 1 dp D^ ^ 1 ^ Dw _ 2 1 dp 

Dt'^ p Sx ’ Dt p Dt p Sz^ 


where 


D 

Dt 


d f) S 

-- + u ^ +v- 
at ax ay 


+ w 




u, V, w being the components of velocity and p the pressure at 
{x, y, z). The equation of continuity was also given by him. 
It was first published in the Memoirs of the Academy of Sciences 
of St. Petersburg in 1770. It is 


dp d(pu) d(pv) d(pw) __ 
^ Sx dy Sz 
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With these equations Euler proved Torricelli’s theorem in the 
following manner : — 

When the motion is steady, we have — = 0. Also, 

because the only external force is gravity, we have X=F = 0, 
and The equations of motion become 


du Bu Bu I Bp 

^ T- + ^ + ?/; tt” = t 

Bx By Bz p Bx 

Bw Bw 
Bx By 


Bv Bv Bv I Bp 

Bx By Bz p By 


Bw 


1 ^ 
p Bz 


Since p is constant for ordinary fluids, the equation of con- 

. . , Bu Bv Bw 

tinuity becomes ~ ~ = 0. 

Bx By Bz 

Further, since u, v, tv are the component velocities in an irro- 
tational motion, we have 

BS B<f> B<f) 

Hence ^ ^ . And we have two similar equations. 

By ByBx Bx 

Bu Bv Bw A 1 ^ A Bp 

“s*" s 


SimiUrly, J ^ .„d j i 


p Bx 

I Bp ^ I ^ A Bp 

- ^,and J (q^) =g 

p By Bz p 


Multiply these equations by dx, dy, dz respectively and inte- 
grate; then 


which is a relation between the pressure and velocity. 

Since the hole is small, we can neglect g* at the level of the 
water. This gives 


0 = 




? + C’. 


.-. i g* = —11 — ? + gz. But ^-11 — ^ = 0, at the hole ; 

P P 

at the hole, J q^-=^gzy which proves the theorem. 

• It will be noticed that the equation of continuity, viz., 

^ f ^ ^ = 0 becomes, for irrotational motion, ^ ^ ^ ^ 
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which is the same as Laplace’s equation about which I spoke 
yesterday. Problems in Hydrodynamics can therefore be treated 
in the same way as problems in Electrostatics or Heat. 

For rotational or vortex motion, the equations of motion 
remain the same as before. We have also 

where t?, £ are the components of the curl. 


Then 


Dt\p} p dx p Sy^ p dz ' 

D/r)\^^dv 7) Sv ^ dv 

Dt\p) pdv^pdy^pdz^ 

Dt\p) pdx^pSy^pSz' 


We can solve these equations in terms of a scalar potential 0 
(no longer a velocity potential). 


If 0= ” ^ dx' dy* dz\ ~ dx' dy' dz\ 


etc., 


r being 

and 

then 


s/{x^x'f ¥{y--yy^{z^zy . 

dv dw 
~ dx^ dy^ dz * 

- ?? 

dx ^ dy dz ’ 

dy ^ dz dx ' 

dG^dJl 
^ dz ^ dx dy ' 


We can thus express u, v, w in terms of F, 0,H. The anal- 
ysis is very much the same as that which is used for problems on 
the motion of electricity. This is an important point of view 
and is of special importance in the treatment of viscous fluids. 

The first attempt to study the motion of viscous fluids was 
made by Navier in 1827 and the investigation was extended by 
Poisson in 1831. But all the investigations which treated the 
problem from the point of view of molecular theory, took the mole- 
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cules as mere points, although they were supposed to be of finite 
size at the start. These investigations were therefore only fruitful 
in creating interest, but Iweie not of much scientific importance. 
The first person to dea with the problem with any rigour was 
De St. Venant, who is famous for his work in the theory of 
Elasticity. The equations of motion are 


Du \ V dB 

Dt p dx^ % dx^ 





I dp V flB 
p dy 3 dy 


+ V r ‘V, 


Div \ dp V dB , 

W~ ~p 

where Q is the dilatation and equals ~ ^ ^ . 

^ dx dy dz 


These equations have a meaning only when the first and the 
second differential coefficients of ii, v, w exist and arc finite. The 
continuity of these had also been postulated, in deriving the equa- 
tions. Prof, Oseen showed that the motion can be determined even 

in the case when ^ ^ , etc., do not exist at all. He 

dx^ dt' dx^ 

critically examined the equations, just as I have done for the 
conduction of heat * and obtained 

j { nx, t) - T(x, 0) ]dr = ^ 

-IT 0 


as the equation of the phenomenon of linear conduction. 

Oseen has obtained certain very interesting results. I read 
out to you a translation of his own views on the subject : ‘‘The 
differential equations of Navier admit of being transformed into 
integral equations in which the differential coefficients 

du dv dw d^u d^u dp dp dp 

W ’ di. ’'m ’ ^ 'df ' ’ 

do not occur. Indeed it can be shown that the integral equa- 
tions possess systems of solutions for which these differential 
coefficients do not at all exist. One must then ask whether these 


t “ Constitution of Matter and Analytical Theories of Heat” {Ahhandbungen 
der Oeaelhchaft der W iasenachaften zu Gottingen, Bd. 2, 1903). 

2 “ Ueber die Bedeutung der Integralgleioh ungen in der Theorie der Bewe- 
gung einer reibenden, unzusammendruckbaren Flussigkeit ” {Arkiv for Matematik, 
Aatronomi och Fyaik, Bd. 6, 1910). UoL' dCZ 
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solutions of the integral equations, which are not at the same time 
solutions of the differential equations, possess only a mathematical 
importance or whether some among them have a physical import- 
ance. I undertake the task to show that the latter is the case and 
that therefore not the differential equations of Navier but the 
integral equations should be considered as expressing adequately 
the physical hypotheses which lie at the basis of the theory of the 
motion of a viscous fluid. 

“ The direct derivation of the integral equations of Hydro- 
dynamics without the intermediation of Navier’s differential equa- 
tions, appears to me to be of interest from two points of view. 
First, obviously the theory gains thereby in generality, secondly, 
the theory gains much in simplicity. In order to see this, we need 
only consider how complicated the considerations are which will be 
necessary to decide under what circumstances the potential func- 
tions will possess differential coefficients of the second order. 
These questions and others of the same kind occur in Hydro- 
dynamics if one starts from Navier’s differential equations. If one, 
however, obtains directly the integral equations of Hydrodynamics 
all these questions lose their importance. A similar simplifica- 
tion is obtainable in a similar manner in some branches of mathe- 
matical Physios. I may mention as examples the theory of poten- 
tial, and the theory of the conduction of Heat.” 

I give you now the integral equations of Oseen. They bear a 
resemblance to the equations for the vibrating string with an 
infinite number of edges. 

From the principle : 

Amount of motion coming into a closed volume 

= amount of motion due to extraneous forces, 

we get 

P H dr^X'ir) dv + dryX. (t) ds = 0 ; 

(O ^ (M) f S 

and two similar equations ; 

u V/ / du du Su\ 

where fa — + 

V dr dy / 

and X„ = Xg, cos nx ^ Xy cos ny + X^, cos nz. 

Here A*, Xy X, are stresses acting in a plane perpendicular 
to the ir-axis. In the case of a perfect fluid the stress is parallel 
to the axis. In that case Xy and A. are zero. Thus in the case 
of viscous fluid there are nine component stresses to consider. 

The differential equation has thus been replaced by an inte- 
gral equation. There is no great difficulty in the solution of such 
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an equation and we can show that the old solution satisfies this 
equation. 

T now come to Electrodynamics. The equations are 


Sy fijS 



dx ’ 


4t7TW = 


^__r9a 
dx dy* 


dt dy dz 


dt ~ dx ’ 


dt'' dx dy ’ 

where u, v, w are the components of the total current strength at 
(a;, y, 2 ), a, ft, y are the components of the magnetic force, a, &, c 
are the components of magnetic induction and X, 7, Z are the 
components of the electric force. The equations can be solved 
by the use of vectors (F, 6?, H) and the potential 0 (cf. rotational 
motion in Hydrodynamics). 

The equations of the theory of electrons arc more general, 
and include the equations given above as a particular case. 
Prof. H. A. Lorentz has studied the motion of electrons. The 
equations of motion are 

IdX u_dy dj 
c dt c~ dy dz ' 

1 dY V da dy 
c dt ^ ^ dz dx * 

IdZ w d ft da 
c dt ^ ^ c dx dy ’ 

dx ^ dy^ dz 

da dZ dY ^ ^ , 

where p is called the volume density of the electricity and u, v, w 
are regarded as the component velocities of the element of electricity 
which is at the point {x, y, z) at time i. 
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The method 
Put 


of solution is 

as follows 

a 


BG 


^ By 

-Bz 

t 

p 

_^BF 

^BH 


Bz 

Bx 

* 


BG 

BF 


y 

^Bx 

^Sy 


Y 

IBF 

BV 

A 

( 

; Bt 

Bx ’ 

V 


1 BG 

BV 

I 


c Bt 

"By ' 

r/ 


IBH 

BV 

/j 


c Bt " 

Bz 


where, F, G, //, are solutions of 


1 iVF u ^ 
r -TT-r +P ~ = 

dl^ ^ c 

1 V 


c' dl^ 


and 


c 

w 


vv/- i/^ + p-=0; 

c* dr c 


also 

and 


2Tr 1 

^ ^ ^ 1 

^ By ^ dz ^ V. dt 


There is no difficulty now iu getting a solution. If p-=0, we 
have the equation of wave propagation. From this solution as a 
basis we can build up the complete solution. 

Partial differential equations are found to bo quite serviceable 
for all ordinary purposes. But in a rigourous treatment of any 
problem they have to be relegated to a secondary place. The 
work of the last 120 years is of importance as it shows that 
most results and processes dealing with differential coefficients are 
valid only under certain conditions. Thus Taylor’s theorem is not 

1 

always valid. For example, the expansion of e according 

to this theorem, viz.. 

i 
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is wrong because every member of the expansion is zero. This 
discovery, due to Cauchy in 1815, astounded Lagrange who believed 
that if the differential coefficients of all orders existed and were 
finite, the function would be expansible by Taylor’s theorem. But 
I have mentioned to you an example which is in direct contradic- 
tion to his belief. However, in spite of its limitations, the partial 
differential equation will continue to be useful. 

Now I will briefly recapitulate the essence of what I have said 
in criticizing the common use of partial differential equations. 
If we have a string with an edge, the first and second differential 
coefficients do not exist at the edge. The equation of motion be- 
comes meaningless in this case, [f we have the initial state of heat 

f{x) = r for x>{), and f{x) - -x f or .r < 0, 


the first and the second differential coefficients are non-existent at 
the origin. Consequently, the equation of the linear conduction 


of Heat 


has no meaning for t = 0, x=.0. 


In Electrostatics, 


Magnetostatics or Gravitational Potential, if /3=:cos log | 

does not exist. In Electrodynamics too a similar remark applies. 
In Hydrodynamics, as Prof. Oseen says, it is not always the 
differential equations but the integral ccpiations which truly express 
the fundamental fdiysical hypotheses. 

It is quite possible that at some time in the future we shall be 
so advanced in our thoughts that a diffei'ential equation will appear 
to us as but a very crude mathematical instrument, and we shall 
discard it in favoui of the more powerful and more refined integral 
equation . 

It has been a real pleasure to me, gentlemen, to come and 
deliver lectures before you, and 1 shall be much pleased if I find 
that luy lectures have been of use in fostering a spirit of research 
here. 


Baptist Mission Press, Calcutta. 








